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1  .  Introduction' 

In  [223,  Plato  states  that  his  teacher,  Theodorus  of  Cyrene, 
was  th<^  first  individual  to  prove  that  -/3,V5,.  .  .  ,  Vl  7  were 
irrational  12,  p. UbJ  .  Since  then  many  scholars  have  wondered 
how  he  did  this  with  the  mathematical  tools  he  had  at  hand, 
and,  perhaps  more  puzzling,  why  he  stopped  at  the  number 
seventeen.  An  answer  which  is  almost  certainly  incorrect  but 
which  is  nonetheless  intriguing  was  given  by  J.  H.  Anderhub  in 
1918  Csee  [3,5,6,10] J.  He  suggested  that  Theodorus  might  have 
used  the  following  construction.  Form  the  triangle  with 
vertices  at  CO, 03,  CI, 03,  CI,  13.  Note  that  the  hypotenuse  has 
length  -/Z.  Starting  now  at  the  point  CI, 13,   draw  a  line 

segment  of  unit  length  perpendicular  to  the  hypotenuse  and  in 
the  counterclockwise  direction;  this  gives  a  new  triangle  with 
this  line  as  one  side,  the  previous  hypotenuse  as  the  other 
side,  and  a  hypotenuse  of  length  V3.  Continuing  in  this  way, 
a  spiral  figure  Ccalled  the  guadratuwrzelschnecke  [103,  or 
square-root-snail  3  emerges  with  the  numbers  Vn  as  the  lengths 
of  the  radii  CFig.  13.  The  last  hypotenuse  drawn  before  the 
figure  overlaps  itself  has  length  Vl7,  and  this,  it  was 
suggested,  was  why  Theodorus  stopped  at  seventeen  Csee  C5D3. 

Although  this  is  almost  certainly  the  wrong  answer,  the 
sequence  of  points  formed  by  this  process  was  studied  on  its 
own  merits  by  Hlawka  CC10];  see  also  [£4]3.  He  proved  a 
number  cf  results  concerning  this  construction.  The  study  was 
then  picked  up  by  Davis  C3],  who  noted  that  the  points  could 


2  =     Z         +     iZ     /|z       I  CI. 15 

ri  +  i  n  on 

Cwith  z  =1  ~)  and  that  they  were  asymptotically  an  Archimedean 
spiral.  He     then     interpolated      an     analytic      curve     to     these 

points,     which    he    called    the    Spiral    of    Theodorus. 

Noting    that    CI. ID    could    be    considered    Euler's    method    for    the 
ordinary   differential    equation 

z    =    z/|z | 
with    unit    step    size,     Davis    suggested    as    a    generalization    of 
CI. ID    the    iteration 

z  =    otz      +    /?z   /  I  z    I 

r>+i  n  n  n 

with     a>fi^<£     and     z        a     given     complex     number.  This     iteration, 

called  the  Complex  Generalized  Theodorus  Iteration  CCGTTD,  is 
studied      in      [3,15,201.  Davis      later      suggested      the      further 

general i zati  on 

V  =    A*V      +    B*V  /IIV    II  CI.  S3 

where  A  and  B  are  real  mxm  matrices,  V  is  a  qi  ven  nonzero 
m- vector  ,  and  II  •  I!  is  the  Euclidean  vector  norm.  This 
iteration,  named  the  Generalized  Theodorus  Iteration  CCliJ  by 
the  present  author,  displays  a  great  variety  of  strange 
attractors  for  appropriate  choices  of  A  and  B  CFigs.  E-5J  ,  and 
includes  the  complex  GTI  as  a  special  case;  it  is  studied  in 
[4,15,18,191.  In      this      report      will      summarize     some     of      the 

information  contained  in  [181  and  [191  Cadding  some  material 
from  [151D  with  respect  to  sufficient  boundedness  conditions 
and     asymptotic     estimates     for      this     iteration.  Since     we     are 

interes    ed     in     boundedness      and      asymptotics,       in      stating      our 


result?  we  will  ignore  the  case  where  V   is  such  that  V   is 

o  n 

the      null       vector      for      some      n      Cin      which      case      the      iteration 
stops} . 


2.       BOUNDEDNESS 

Let     us     try    to    find    sufficient     conditions    for     boundedness    of 
the    GTT  .        Taking    norms    on    both    sides    of    CI. 2D    gives 

II V        II    <    HAH  •  II V    H     +     IIBH 

n-t-l  n 

Iterating    this    relation,    we    find    that 

n-  l 
II V    II    <    II  All"  •  II V    II     +    IIBII-r    IIAII*  C2.1D 

n  O 

t  =0 

and    clearly   if    II A  IK  1    then    in    the   limit    the   norm   is    bounded    by 

IIBH  /CI    -    II  All  D 
Hence     I'AIKI     is     sufficient     for     boundedness     of     the    orbit     Cfor 

any    V    D.         From    this    and    the    similarity    of    C 1  .  2D     to    the    linear 


■ysiem 


V  =    A*V  C  P. .  2D 

n-t-l  n 


which  has  a  bounded  solution  for  every  V   when  the  spectral 

o 

radius  of  A  is  less  than  unity  Cand  at  least  one  unbounded 
solution  when  the  spectral  radius  of  A  is  greater  than  uni tyD , 
it  seems  reasonable  to  conjecture  that  pCAD<l  suffices  for  the 
boundedness  of  C 1 . 2D ,  where  pC AD  is  the  spectral  radius  of  A. 
That  this  is  in  fact  the  case  follows  from  two  lemmas.  We 
omit  the  rather  technical  proof  of  Lemma  1.  It  is  not  hard  to 
establish  this  result  by  using  an  approach  similar  to  that 

TteoH         i    r-i        c-car-4    -t   r-M~>^-        f  i   *\/-c.        onH        ^-  i   *v        Vm  t4  i    +  t   c         rather         tedi  OUS        tO 

carry     cut     the     proof      in     detail      Csee     C153D.  Lemma     2     is     a 

discrete    version    of     an    ordinary    differential     equations    result 


f ^.tf-iH     i  *-,     r  1   1 

Lemma    1  :        If     II V    II     is    unbounded    -then     II V    II — too. 

n  r> 

Lemma    2  Let     A    be    a    matrix    with    pC AD < 1     and    suppose    that    BC  n} 

is    a    sequence    of    matrices    with     HBCnDll — >0    as    n — >oo.         Then    all 

solutions    of 

z  =    CA    +    BCnDJ  -z 

tend    to    zero    as    n — >oo. 

Proof:        Since    pCADvl     there    exists    a    matrix   norm    II  •  II       such    that 

A 

II A II    <1.        We    will     choose    this    to    be    a    natural     matrix    norm    and 

A 

will     also    denote     by     II  ■  II        the     vector     norm     which     induces     it. 

A 

Let      e      equal       CI -II A II    ~)  .  Now      since      all       matrix      norms       are 

A 

equivalent    in    the    sense    that    if     IIMCnDII    — >0    then    II MC  n~>  II  _ — >0,     we 

o  ft 

have    from   our    hypothesis    that     IIBCn3H    — »0.        Hence    we    can    choose 

A 

an    N    such    that     llBC  O  II    <e/S    for    all     i>H. 

A 


z  -    C  A    +    BCnJ  'J  -z 

n+i  n 

Cfor    sorre    z    D     N    times    to    produce    the    vector    z  Set    to    =z       and 

O  N  O  N 

(•<■  =    CA    +    BCN    +    n))  -w 

r>+l  n 

Vn>0.        Clearly   to  =z  for    all     i>0.        Now 

l  N-K 

go  =    CA    +    BCN    +    nlO*     ■  •  ■     *C  A    +    BC  N30  *co 

n+i  O 

and    so 

llw        II       <    HA    +    BCN    +    nD  II    *    ■•■     *IIA    +    BCNDll    *llw    II 

n-t-1     A  A  A  O     A 


Then 


CO  II 

nt-i     A 


O      A      '  ' 
v  =0 

A 

n 

<  ii  <*>  ii     n  c  ii  a  ii 

O     A      ■  '                     A 

x.  =0 
n 

+     II BC  N    +    i  5  II     Z> 

A 

<  iiw  ii       n  cci 

O     A          '• 

v  =0 

r> 

-  <=:>   +  <=/2:> 

=  ii  co  ii       n  ci   - 

O     A          '  ' 
v  =0 

^y"ED 

whi  r.h    clearly    lends    to    zero    as    n — too.         Hence     llco    I!    — >0    as    n — >a" 

n     A 

so    that     ilz    II    — >0    as    n — >oo,     and    so    by   the    equivalence    of    vector 

n     A 

norms    on    [R       we    have    that     Ilz    II — >0    as    n — >oo.         This    completes    the 

pr  oof .     d 

We    may   now   show   the   desired    result. 

Theorem    1:        If    pC  A3  <  1     then    the    solution    of    CI. 2D     is    bounded. 

Proof:        Suppose    the   contrary   and    choose    a    V      such    that     II V   II    is 

O  n 

unbounded.    By    Lemma  1,  II V  II— >oo.    Define  BCn5=B/HV  II.    By  the 

n  r> 

hypothesis,     HBCnDll — ►O.        Consider    the    iteration 

z  =    CA    +    BCnlO^z 

n-M  n 

with    z    =V    .        Clearly   z    =V      for    all    n.        The    conditions    of    Lemma 
o       o  o        o 

2    are     met     and     so     II V    II — >0    as     n — ><x>.         But     this     contradicts      IIV    II 

n  n 

unbounded,     which    establishes    the    result,     d 

It    i^    also   possible    to    show   that    if    pC  A')  >  1    then    there    exists 
a     V       such     that      HV    II     is     unbounded     C  see     section     5D  .  In     this 

O  n 

sense     the     GTI      is     surprisingly     similar      to     the     linear      system 

C  cZ.  C.J. 

In    the    particular    case    A=0    a     detailed    analysis    is     possible 
[16,37],  This     case     is     related     to    the     power     method     for     the 

numerical    solution    of    eigenvalue    problems. 

3.       ASYMPTOTICS 


This  leaves  the  question  of  boundedness  for  the  case  pC  A~J  =1  . 
It  is  possible  to  construct  examples  showing  that  when  A  is 
uni  spectral  the  orbits  of  the  GTI  may  be  always  bounded  for 
one  choice  of  B  and  always  unbounded  for  another  choice  of  B. 
However,   a   general   statement   about   the  asymptoti cs   of   the 


F, 


4  4  m^mJ]   .  f-',^,     Wat*    ViCi     fffaHca    fhci?--:     A     -j  >-     ^.4"     hnttfli^^H     t  ^•"M^        i  H^  +      1  *^        'i'hpf; 

pCAD<l      and     any     eigenvalues     of     A     'that-     lie     on     the     unit     circle 

are         simple         Ci.e.  have         as         many         linearly         independent 

eigenvectors    as    their    mul tiplici tyj .        In    this    case    the    linear 

system    C2.2D     has    the    property    that    all     solutions    are    bounded 

C although     in     general      the     bound     will      depend     on     V    D.  Before 

stating    the   corresponding   result,    for    the   GTI  ,    we   state: 

Hukuvara's    Theorem?       Let    z     be    any   solution    of 

t 

t       =  ca  +  ect;o*z  C3.  id 

for    t=0,l ,2,. . .     where    A    is    a    constant    matrix   and 

£    HBCtJ  II    <    co 
t  -  i 

Suppose    that     A    is     of     bounded     type,     so    that     all     solutions     of 

y  =A*y  are     bounded     for     t>0.         Then    every    solution    of     C3.1D 

is    bounded. 

Thi  s  result,   an  analogue  of   a  theorem  from  differential 

equations,  mav  be  found  in  [213.    The  proof  there  uses  the  I 

l 

vector    norm,     but    the    argument    is    valid    for    any    vector     norm    and 
compatible      matrix     norm.  We      will       now      use      it      to      prove      an 

asymptotic    result    for    the   GTI. 
Theorem    2:        Let    A    be    of    bounded    type    in    C 1 .  2D  .     Then 

00 

V   1/IIV    II     =oo 

Proof:        Suppose    the    contrary,    i.e.     that    for    some    V      we    have 

oo 
T   1/IIV    II    <    co  C3.  2D 

n 
n  =  0 

and    construct    such    a    sequence    <V  >    by   using    this    V      and    CI. 2D. 
Define    BCnD=B/IIV    II    and    consider     the    iteration 

Y  =    CA    +    BCnDD*Y 

n+l  n 

where      ve     set      Y    =V    .  Clearly     Y    =V       for      all      n>0;      the      two 

O         O  n         n 


sequences    &v e    identical.        Now 

E    HBCrDII     =       £    HBII/IIV    II 

n 
n=0  n=0 

00 

=    IIBII    £  1/IIV    II 

n 
n  =  0 

which      by     supposition     is      bounded.  Therefore      by     Hukuwara's 

Theorem,     II Y    II    is    bounded.        Then    II V    II    is    bounded,     that    is,     3M>  0 
such     that      II V    IK  M     for     all      n.  But     then      111  /V    ll>l/M     and     since 

n  n 

1/IIV   II— x-»0       we       cannot       have       C3.ED.  This       contradiction 

establishes    the   result.     □ 

iii  fact  a  slightly  more  general  result  holds.  If  II*  II  is  any 
vector  norm  and  K'J  is  any  function  that  is  continuous  on 
rO.ocO  and  positive  on  C0,oc0,  and  if  A  is  of  bounded  type,  then 
any   solution    of 

V  =    A*V      +    B*V  /$C  II V    ID 

n+l  r>  n  n 

must    be    such    that 

on 

£   1/$C  II V    IO     =oo 

n 

Thi  s     may    be    shown    in    a    completely    analogous     manner     Csee     [15] 
for    details    and    applications!). 

Theorem  £  establishes  that,  when  pC A3 =1  but  A  i s  of  bounded 
type,  if  the  orbits  diverge  then  they  do  so  in  a  way  that  is 
not.  much  worse  than  linear  C  si  nee  the  sum  of  the  reciprocals 
of  the  normc  rii  verges).  In  fact  if  HAll=l  in  the  spectral  norm 
then  it  is  a  fact  that  the  divergence  is  at  most  linear,  as 
can    be    seen    by    setting    llAll=l     in    C2.13    to    get 

II V    II    <    II V    II    +    n IIBII 
as    the    asymptotic    estimate.        We    will    later    consider    this    point 
in    greater     detail     when    we    show    in    section    5    that    if     A    is    of 
bounded    type    then    the    divergence    is    always    at    most    linear. 


4.    Green's  Function  Representation 

In    order     to    refine    our     asymptotic     estimates    of     the    previous 

section,       we      will      need      a      formula      for      the      solution      of      the 

iteration    C 1 . 2D .         Given    an    iteration    of     the    form 

Z       =    ACtJ*Z  +    w 

t  i-1  t-i 

with    o.%      and    z       given    vectors     and    ACtD     a    time-varying    matrix, 

the    solution    can    always    be    written    in    the    form 

t 
z    =  r  Y  -Y^-w  +  Y   -Y"*-c 

t  +  1  t    S      S        t  +  1     o 

S=0 

where  c  is  the  initial  condition  z  ,  and  CY  >  is  a  fundamental 

o  t 

matrix    set    of     solutions    for 

z      =    AC  1 3*7  C4.  ID 

t  t-i 

Cthat   is,   Y  =1   and  Y   =ACnD*Y  D   CS13.     For   the   GTT  ,   the 

O  n+l  n 

matrix    A    is    constant,     and    z    =V       is    given.         Furthermore,     Y    =A 

o        o  n 

for    the    linear    system    C£.£D    corresponding    to    C4.1D.,     and 

w      =    b*V  /IIV  il 
t  t  t 

so    that 

n 

V  =    An+1V       +       r   An_t*B*V/iiV   II  C4.2D 

n+l  O  "*  t  t 

t  =o 

is  the  solution  of  CI.  ID.  Thi s  is  called  the  Green's  Function 
Representation  of  the  solution.  We  will  use  this  form  of  the 
GTI    to    address    the    questions    of    boundedness    and    asymptotics. 

5.     A    Direct    Approach    to    Boundedness 

The  approach  used  in  the  second  section  to  show  that  pCAD<l 
suffices     for     the     boundedness     of     CI. 2D     is     the     one     originally 


employed     by    the    author     in     [15]     and    also    appears     in     [183.         A 

more    direct     method,     however,     was     used    by    the    author     in     [193, 

and    we    now    give    this    method. 

Suppose    that     pCA3<l.        Then    there    exists    some    natural     matrix 

norm     II  •  II        such     that      II A II        <     1.  Using     II  ■  II        to     represent     the 

ex  a  a 

vector     norm     which     induces     thi  s     msi.rix     norm     as     well,      we     have 
from   C4.23    that 


II V        II       <    IIAn+1V    II 

r.+  l     CX  O     CX 


n 

r  An-t*B*V  /iiv  I 
*"  t         t 

t  =o 


ex 


Vn>0.        Thus 


n 

II V        II       <    IIAIIn+1-HV    II       +       T    IIAI^-IIBII     -k 
n+i    c*  ex  o   a  a  a 

t  =o 

where  >  is  a  positive  constant  such  that 

II VII    /-IIVII    <    k 

CX 

for  all  nonzero  vectors  V.  Such  a  k  exists  since  all  vector 
norms  are  equivalent  on  D?  [123.  Summing  the  geometric  series 
and    collecting    like    terms    gives 

II V       II       <   eCrO    +    k'itBll    /CI    -    II A II    3 

n+l     C  CX  O 

where    c-.  n.l'     >0    as    n — »ao.        Hence    II V    II       is    asymptotically    bounded, 

n     CX 

and     bv      the      equivalence      of      vector      norms      so      is      II V    II.  An 

n 

asymptotic    bound    for     II V    II       is    given    by 

n     Cx 

k-UBll     /CI     -    II  All    3 
ex  a 

in  terms  of  the  cx-norm.  When  the  spectral  norm  of  A  is  less 
than    unity   the    corresponding    bound    using    that    norm, 

IIBH/C1     -    IIAII3, 
is    generally   fairly   tight. 

We  have  shown  that  pCAIKl  suffices  for  boundedness.  In  an 
anal  ago. is  way  we  can  show  that  if  pCA3>l  then  there  must  exist 
some    V      such    that     II V    II     is    unbounded.         CA    similar     approach    has 


10 


hfcjcifi      ti   -fil  Mtfea'J      1--.*      V    — t 


Chiton :     <pa^    r  1  ^i'i 


For,     considering 


C4.  S3  ,     we    have    that 


V         I!     > 
n+1 


A^V 


r  A       *B*V  /IIV 
^  t        t 

t  -o 


C using    ;.he    spectral     nornO.         Let    X    be    an    eigenvalue    of     A    such 

that       \\\=pQAD       and      choose      V        to      be      an      eigenvector       of       A 

o 


associated    with    the    eigenvalue    X.         Then    we    h 


V        II     > 

n+1 


\mlV 


ave 


T  An  l*B*V  /IIV 
*"  t  t 

t  =o 


and    since    the    norm    of    the    sum    is    at    most    the    sum   of    the    norms, 


V        II    > 

n+1 


|x  r+1nv    II    -      r    II  A  II l-  II  B  I 
o 

t  =o 

p     |X|n+1    -     IIBH-C1     -     HAlin+1D/Cl     -     IIAIO 
o 


wi  t  h    p    =  li  V    li 
o  o 


then    c>0    and 


I f    we    1 et 

c    =    IIBII/C  HAH    -    ID 


V       II 

n+1 


p      X  -    c  II  All  +    c 

o 


in    term?    of    c.        Now    pC  AJ  =  |  X  | -t-e    for    some    e>0,     so 
II V        II    > 


n  +  1 


p     I  X |  -    c  C  |  X  |     +    eD  +    c 


X|n+1Cp       -    cCl     +    €/|Xpn+1    +c 
o 


C5.  13 


If  e  n;  equal  Lo  zero,  we  may  now  choose  p   to  be  or eater  than 

o 

c    in    order     to    have     II V    II — >co.         Otherwise,     the    term    involving    e 
will    grow   arbitrarily   large,    again    giving    II V    II — ►co. 

n 

Again,       we      remark       that       the       GTI       Cl.iJ       has       boundedness 

properties     similar     to    those    of     the     linear     system     C2.23.  In 

the    next     section    we     will     derive     improved    asymptotics     for     the 

indeterminate    case    where    A    has    unit    spectral     radius.         Before 

doing    so    let    us    collect    the    information    above    as    a    theorem. 

Theorem    3:         Consider     the    iteration    C 1  .  2D  .         If     pC  AD < 1     then     the 

iteration     is     bounded     for     any     V    ,      and     if     pCA3>l     then     there 

o 
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exists    «    V      bocii    chat    the    iteration    is    unbounded, 
o 


5.     Improved    Asymptotic    Estimates 

From  C5.13  it  is  immediately  apparent  that  when  pCAD>l  the 
divergence      will,       in      general,       be      exponential.  This      again 

leaves  the  case  pCAJ=l.  As  before,  let  us  first  look  at  the 
case  where  A  is  of  bounded  type.  Using  the  spectral  matrix 
norm,     C4.83)    gives 


V        li     <     IIAn+1V    II     + 

r>+l  O 


r. 


'   A^  l*B*V  /IIV 
**  t  t 

t  =o 


C6.  ID 


<    K      +       T    H  A*  II  •  II  B  II 
1  *"* 

t  =o 

where      K  >0      is      a      bound      on       II A        V    II       C  which      in      general       will 

depend    on    V   }  .         Now    since    A    is    of     bounded    type,      II A    II     is     also 
o 

boundeo.     that    is,     II A    II  <K    for    some    K>0.        Thus 

n 
II V       II    <   K      +    IIBH    r  K 

n+1  5 

t  =o 

<    K      +    K  n 

1  2 

showing    that    the    divergence    is    no    worse    than    linear. 

If     pCAD=l     but    A    is    not    of    bounded    type    then    it    is    possible 

to    show    that     II  A^li  =<DC  n™}     Cin    fact     II  Anll  =©C  np-1D     where    p<m    is     the 

maximal     degree    of    all     nonlinear    divisors    of     A    associated    with 

unimodular    eigenvalues!),     where    A    is    mxm.        Hence    from    C6.13     we 

have 

n 
II V        II    <    II  A^"1  II  •  II V    II     +       F    II  A1  II-  IIBII 

n+l  O  *"* 

t  =o 

and    clearly    II V       II     diverges    at    most    as    OCn       D.         Cin    fact,     it 
diverge-     as     OC  n   D  ,      and     hence     in     the     worst     case     as     £>Cn    D  D  . 

Hence,     when    pC  A"J  =1  ,     if    the    iteration    diverges    then    it    does    so 


IS 


O.  •••         lil'-',---   ■-  l^'.-ij.    ^  i  I'.'lHi    O.J.    J.   J    ■  "'■  *»  '      •*■     '-•'-•  S.-4  1J.    w         CA  -5         O.  ^il'.'^l     »_•  III. 

Theorem     4:           Consider      C 1 . £j      and     suppose     that  pC AJ =1      and     the 

iteration     is     unbounded.           If      all      eigenvalues  of      A     with     unit 

modulus      are     simple     then      the     divergence     is  no      worse      than 

linear.        Otherwise    the    divergence    is    no    worse  than    polynomial 
of     order     m,     where    mis    the    order     of    A. 

7.     Li  mi  t    Sets 


Let  us  define  I~CA,B,V  ;  to  be  the  limit  set  Cor  c^-l  i  mi  t  set 

o 

[131 D  of  the  iteration;  that  is,  Ye.rCA,B,V  J  if  there  exists 

o 

an    incrR^cinn     subsequence    nCiJ     such     that     V       — >Y     as     i — too    when 

the    initia]     condition    is    V    .        Twn    theorem?    follow    easily    from 

o 

the    material     in    [131. 

Theorem     5:  Suppose     that     for      some     V       the     sequence     -CV  '}     is 

O  n 

bounded  Then      FCA.B.V   D       is      a       non-empty      compact      set.      and 

O  r      .  . 

V  — *rCA,B,V  j     as  n— ►». 

r.  O 

Proof:     Clear]  y  the  closure  of   F  is  simply  T  and  so  T  is 
closed.     The   boundedness   of   <V  >   implies   that   T   is   also 

n 

bounded.  Hence      T~       is       compact        C  by      definitions.  By       the 

Bol zanc-Wei erstrass    Theorem    [£33,     T    is    not    empty. 

Now      v  — >r      means      that      inf  Cli  V  -Y  II  t  Yerj  — >0      as       n— *a>      CI  33. 

ri  n 

Suppose  the  contrary.  Then  we  can  find  a  subsequence  of  -CV  >, 
indexed  by,  say,  nC i D ,  which  remains  a  finite  distance  from  T. 
Now      V  is      a       bounded       infinite      sequence      and      so      by      the 

Bol zano-Wei erstrass     Theorem    it     contains     a     limit     point.  But 

that    point    must,    be    in    T,     which    is    a    contradiction,     d 
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Theorem  6:     Suppose  that  for   some  V   the  sequence  <V  >  is 

O  ^  n 

bounder!       and       OerFCA.B.V   "5.  Then       FCA.B.V   J        is       positively 

o  o  r  J 

i  nvari  ant  . 
Proof:        Define 

TCV}     =    A*V    +    B*V/IIVII  C7.1D 

so    that 

V  =    TCV  D 

is       the       iteration       under       consideration.  Since       0«2T,        the 

transf crmati on     T     is     continuous     in     some     neighborhood     of     any 
point     VeF.  Fix     Y     and     choose     a     subsequence     nCiD      such     that 

V       — »Y.        By    the    continuity    of    T,     we    have    that 

tcv     :>_>.tcy:> 

as    i  — ►<*      But 

TCV       j     =    V 

r><v.>  n<i>+l 

so    that 

V       — >TCYD 

ttKv.  > 

where       mCi D =nCi D +1  .  Hence       TCYDeP,        sod       r       is       positively 

invariant,     d 

I  n    f  a  ■  t ,     .i  i 

and  TCYD  is  defined.  In  other  words,  F\<0>  is  positively 
invariant.  Further  topological  results  on  discrete  dynamical 
system?  that  can  be  extended  to  this  case  Cwith  some  care 
taken  near  the  originJ  may  be  found  in  [143. 

Suppose   that   YeTCA,B,V  D,    so   that   V   — >Y        for    some 

O  r>a> 

increasing       subsequence       nCiD.  Since       the       transformation 

T:  IF    — ►£       of    C7.1D    has    the    property   that 

TC-VO     =    -TCVJ 
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it     fed  lows     that     the    sequence    generated    by    -V       is     -V        ,      which 

tends     to    -Y ;     that    is,     -YeFCA,B,-V   ~j  .         In    a    similar     wav,     it     is 

o 

easy   to    see    that    -YeTC-A,-B,V  J    and    YeTC-A, -B, -V  3. 

o  o 

In    the    case    of    a    global     attractor,     FCA,B,V  J     is    independent 

o 

of  V   for  the  particular  A  and  B  under  consideration.   In  this 
o 

case,  it  follows  from  the  above  that  the  attractor  is 
centrally  symmetric,  i.e.  YeF  implies  -YeF",  and  that  the  same 
attractor  is  found  when  A  and  B  are  replaced  by  -A  and  -B, 
respectively.  Some  of  the  attractors  displayed  by  this 
iteration  for  various  choices  of  A  and  B  are  shown  in  the 
l  igurtoj  see  also  [4,  153.  I  series  has  found  an  analytical 
expression  for  the  attractors  for  a  certain  set  of  cases  of 
the  form  A=ocB,  a>0  Csee  [41  J.  For  the  case  A=0  the  form  of 
the  attractor  is  also  known  [15,15,173. 


S.  Discussion 


i  rih=' 


Generalized  Theodorus  Iteration  is  in  some  ways  very 
similar  to  the  simple  linear  system  C Theorem  3D  and  yet  it 
contains  some  very  complicated  structures  CFigs.  2-5D  .  It  is 
likely  that  the  attractors  are,  in  general,  strange  but  not 
always  chaotic  Csee  [93  J;  some  computational  evidence  to  this 
effect  has  been  provided  by  James  Heyman  of  the  Naval 
Postgraduate  School  in  his  study  of  the  use  of  chaotic 
discrete  dynamical  systems,  particularly  the  Henon  attractor, 
as  pse  idorandorri  number  generators  Cpersonal   communication!). 


ares.  are'  in  progress.  Due  -to  ■the  present  plethora  of 
definitions  of  the  terms  "chaos"  and  "fractal"  C see  [15]?  it 
is  difficult  to  prove  anything  substantative  about  the  chaotic 
nature  of  the  attractors  but,  as  noted  in  [11] ,  many 
well- known  strange  attractors  are  still  awaiting  analytical 
verification  of  the  observed  qualitative  features  of  the 
underlying  iterations  and  differential  systems.  One  aspect  of 
the  GTT  which  we  feel  is  particularly  relevant  to  the  field  is 
the  vast  array  of  geometrically  distinct  attractors  to  be 
found  in  the  GTI  for  different  choices  of  the  matrices  A  and 
E;  even  though  these  two  matrices  essentially  represent  eight 
parameters-    the   number   of   distinct   attractors   is   still 

,-,,_,-.,-  ;  r-  ,-   r-.r-,  =T-.H       1  ■»■-•       -.•-(=       r-.-— ,4         -.,.•-.,-.—.       .-.*"        __  _.+  v.—,-        H -'  »-•— «- •—  ■+  c        r-  ?    -•«-.-«- 

^  i-ii    ^  1    a.  -r  i.  j  ,  L^  ■         01:1         we—        o.  -     fc        J  1  -  ■  l         a  woi    «=r        ^>  j  0.1  io^i  jCi  <a  j.  o>  *w  i    <rr  «_  cr        pi.  at  iaj 

map  with  this  property.  Further  investigation  of  this 
peculiarity  seems  warranted.  Hopefully  the  figures  cited  here 
and  in  [151  will  motivate  research  on  the  question  of  the 
number  of  distinct  attractors  which  can  belong  to  a  single 
iteration,  a--  the  works  by  Gleick  [7],  Devaney  [Sj  ,  and 
Wiggins  [£5,c;5]  have  served  to  help  motivate  and  inform  this 
author  . 
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